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The pulsar synchrotron: coherent radio emission 
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ABSTRACT 

We propose a simple physical picture for the generation of coherent radio emission 
in the axisymmetric pulsar magnetosphere that is quite different from the canonical 
paradigm of radio emission coming from the magnetic polar caps. In this first paper 
we consider only the axisymmetric case of an aligned rotator. Our picture capitalizes 
on an important element of the MHD representation of the magnetosphere, namely 
the separatrix between the corotating closed-line region (the 'dead zone') and the open 
field lines that originate in the polar caps. Along the separatrix flows the return cur- 
rent that corresponds to the main magnctosphcric electric current emanating from the 
polar caps. Across the separatrix, both the toroidal and poloidal components of the 
magnetic field change discontinuously. The poloidal component discontinuity requires 
the presence of a significant annular electric current which has up to now been unac- 
counted for. We estimate the position and thickness of this annular current at the tip 
of the closed line region, and show that it consists of electrons (positrons) corotating 
with Lorentz factors on the order of 10 5 , emitting incoherent synchrotron radiation 
that peaks in the hard X-rays. These particles stay in the region of highest annular cur- 
rent close to the equator for a path-length of the order of one meter. We propose that, 
at wavelengths comparable to that path-length, the particles emit coherent radiation, 
with radiated power proportional to N 2 , where N is the population of particles in the 
above path-length. We calculate the total radio power in this wavelength regime and 
its scaling with pulsar period and stellar magnetic field and show that it is consistent 
with estimates of radio luminosity based on observations. 
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1 INTRODUCTION 

The physical mechanism of pulsar radio emission re- 
mains still poorly understood. Observational evidence (e.g. 
Kramer, Wex & Wielebinski 2000) suggests that it must be 
due to some coherent radiation mechanism. Unfortunately, 
neither the exact mechanism (curvature radiation, free elec- 
tron maser emission, excited plasma waves, etc.), nor its 
origin (inner or outer magnetosphere) are clear. It is never- 
theless widely accepted that radio emission originates within 
a few stellar radii from the stellar surface, and is due to rel- 
ativistic plasma outflowing from the two stellar magnetic 
poles (e.g. Kramer et al. 1997). 

We here would like to propose a new physical picture 
for pulsar radio emission that capitalizes on what we have 
learned about the pulsar magnetosphere when studied under 
the assumption of force-free ideal magnetohydrodynamics. 
The first self-consistent solution of the axisymmetric prob- 
lem (Contopoulos, Kazanas & Fendt 1999, hereafter CKF) 
demonstrated that the pulsar magnetosphere consists of a 
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closed line corotating region (the so called 'dead zone' 
extends up to the 'light cylinder' at distances 
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from the rotation axis (fi* is the angular velocity and 
P* = 2n/Q* the period of the rotating neutron star), and 
two open line regions that originate around the two stellar 
magnetic poles (the so called 'polar caps'). From each polar 
cap outflows the main magnetospheric current, which in the 
case of the aligned rotator may be estimated as 
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Here, £>» ~ 10 13 G and r* ~ 10 6 cm are the stellar polar 
magnetic field and radius respectively (Gruzinov 2004; Con- 
topoulos 2005). The =p sign in eq. depends on whether the 
magnetic axis is aligned or counter-aligned with the axis of 
rotation (our results will heretofore refer to the aligned case 
with the corresponding result for the counter-aligned case in 
parentheses). CKF established that the return electric cur- 
rent flows along the separatrix between open and closed field 
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lines. In fact, they self-consistently determined the distribu- 
tion of electric current as well as the distribution of electric 
space-charge everywhere in the magnetosphere. We would 
like to point out here that magnetohydrodynamics does not 
deal with the issue of how this is materialized with actual 
distributions and flows of charged particles (electrons and 
positrons). This is not really a problem for the largest part 
of the magnetosphere since it is known that the number den- 
sity of the electron-positron plasma that fills the magneto- 
sphere exceeds the above space-charge by a huge multiplicity 
factor (on the order of 10 3 — 10 4 ; e.g. Arons 1983), and all 
that is needed in order to obtain the distributions required 
by CKF is a slight difference in the velocities of electrons 
and positrons. As we will now argue, this may not be true 
in the separatrix between open and closed field lines. 



2 THE PULSAR SYNCHROTRON 

We will discuss here only the axisymmetic case (we plan to 
address the 3D problem in Paper II of this series based on the 
numerical results of Kalapotharakos & Contopoulos 2009). 
Inside the closed line region, B^ — 0, and particles at its tip 
corotate with Lorentz factor 7 equal to 
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where, R w r; c is the (cylindrical) distance from the rotation 
axis. Let us heretofore denote quantities at the inner/outer 
surface of the separatrix with =p indices respectively. The 
toroidal magnetic field B~ = 0, and 
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Uzdensky (2003) and others pointed out that, as the tip of 
the dead zone approaches the light cylinder, the poloidal 
magnetic field B~ diverges as 



Bp 



; 7 



(5) 



E is continuous 
1; see 
11 of 



and Bp = (the Lorenz invariant B 2 
across the separatrix, and E = R£l,B p /c; see Ffig 
also fig. 4 of Kalapotharakos & Contopoulos 2009, fi^ 
Timokhin 2006, and fig. 1 of Spitkovsky 2006). 

Note that particles at the tip of the dead zone corotate 
with the star due to the combined action of the magne- 
tospheric electric and magnetic fields. During that motion, 
they deviate around their corotation radius within distances 
on the order of the gyroradius 
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Here, m e , e are the electron mass and charge respectively, 
and u>b = eBp /"/m e c is the gyration frequency. For a 'cold' 
particle distribution, this distance dictates also the thickness 
of the separatrix layer. We will assume that the closed line 
region extends as close to the light cylinder as is allowed by 



the equations of motiorQ, namely that 
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Eqs. © & © yield 
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V10 13 G7 V33ms, 
One interesting feature of the separatrix is that it is 
negatively (positively) charged with surface charge density 

We can express the number density n of charge carriers in 
the separatrix layer as 
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Obviously, the multiplicity factor k must be on the order of 
unity for the force-free assumption to be valid. Note that 
at the tip of the dead zone, the Goldreich- Julian number 
density ugj is equal to 
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(Goldreich & Julian 1969, CKF), and therefore, 
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We have shown here that the number density of charge carri- 
ers at the tip of the separatrix must be close to the minimum 
number density needed to support the surface charge den- 
sity required by eq. ©. These particles corotate with the 
star at speeds very close to the speed of light. If we now 
apply Ampere's law 

47r dh 
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and if we assume that particles at the tip of the separatrix 
move with velocities /3c, then eqs. (|13p & @ imply that 



& ~ Pp ~ Bp 7 



(14) 



Note that the separatrix is negatively (positively) charged, 
and therefore, ji points along the direction of stellar rotation 
and towards the star. In other words, electrons (positrons) 
in the separatrix follow a spiral motion as they move inwards 
away from the tip of the closed line region, with very small 
pitch angles on the order 0/I/7 (fig. 2). 

CKF (and for that matter all subsequent investigations) 
failed to emphasize the significance of this annular (i.e. az- 
imuthal) electric current. These electrons/positrons emit in- 
coherent synchrotron radiation up to a cutoff frequency 
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1 That is, we will ignore here previous solutions (Contopou- 
los 2005; Timokhin 2006) which suggested that the dead zone 
may end at any distance inside the light cylinder. 
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which corresponds to hard X-rays. We emphasize once again 
that this incoherent synchrotron radiation is due to the 
macroscopic corotating motion of the separatrix electrons 
(positrons) around the light cylinder. We would like to name 
this configuration 'the pulsar synchrotron'. 



3 COHERENT RADIO EMISSION 

What remains to be calculated is the total power radiated 
from such an electron (positron) configuration. The power 
radiated by an individual particle is equal to 

with a frequency distribution (in erg/s/Hz) 

F v = 3.23-^-(^P,) 1/3 (17) 

at low frequencies v <C v c (Jackson 1975). As seen in fig. (2), 
the annular volume where radiating electrons/positrons at- 
tain their maximum Lorentz factor given by eq. (0 has 
thickness and height comparable to d, and contains 

/ R \ 1/2 

N e ~ 2irnr lc d 2 « 8 x 10 25 [j^q ) ( 1S ) 

particles. The total incoherent power emitted by these par- 
ticles is probably too low to be observable. As we will now 
argue, however, coherency effects may become very impor- 
tant in the radio frequency range around 100 MHz. 

It is well known that, in laboratory synchrotrons, parti- 
cles travel in bunches with macroscopic dimensions (~mm) 
along their direction of motion, where the bunch length A 
is determined by the geometrical and operational charac- 
teristics of the synchrotron. It has been observed that, at 
wavelengths comparable to the bunch length (or the length 
of any structure in the bunch), the radiation from multiple 
particles is in phase, giving a radiated power proportional 
to N 2 times the power radiated by an individual particle, 
where N is the bunch population (e.g. Venturin et al. 2005). 
This emission is orders of magnitude stronger than the in- 
coherent power at those wavelengths, which is proportional 
to N. A practical application of this effect is the generation 
of THz coherent synchrotron radiation in synchrotrons with 
X-ray incoherent synchrotron radiation (e.g. Sannibale et 
al. 2003). The pulsar synchrotron is certainly very different 
from an earth synchrotron. Nevertheless, we may also obtain 
a characteristic bunch length A in the pulsar synchrotron as 
follows: separatrix electrons/positrons stay in the annular 
volume of maximum Lorentz factor at the tip of the dead 
zone for a path-length on the order of 

(fig. 2). If indeed particles emit coherently at the above char- 
acteristic wavelength during the time they need to cross the 
region of maximum Lorentz factor, the coherent radiation 
corresponds to about 100 MHz. The annular radiating vol- 
ume consists of 



regions of length A, each one consisting of N = N B /(i par- 
ticles. The total luminosity of the coherent emission in a 
frequency range on the order of c/A may be estimated as 
the sum of the luminosities of the coherently emitting re- 
gions 

L ~ 2fiN 2 vF„l =c/x 

~ 4 x 10 28 erg s" 1 f-^V /6 f-^V 7/3 (21) 
s U0 13 G/ 1 33ms/ v ; 

(the factor of 2 accounts for the contribution to the total 
radio luminosity from both hemispheres of the dead zone 
separatrix) . 

This is just an order of magnitude calculation of the 
coherent radio luminosity expected from the tip of the dead 
zone. More detailed analysis will require the 3D numeri- 
cal self-consistent reconstruction of the electric and mag- 
netic field geometry near the tip of the dead zone together 
with the calculation of separatrix particle orbits (work in 
progress). Nevertheless, our estimate of the coherent radia- 
tion frequency is not far from the observed peak of pulsar 
radio emission at a few hundred MHz. Moreover, as we can 
see in fig. (3), our result in eq. (|21|l is consistent with esti- 
mates of the total pulsar radio luminosity obtained from the 
ATNF pulsar catalogue^ (Manchester et al. 2005). We there- 
fore believe that our scenario merits serious consideration as 
an alternative to the radio pulsar phenomenon. 

4 DISCUSSION 

The scenario that pulsar radio emission may be coming from 
the light cylinder is not new. In the early days of pulsar as- 
tronomy, people discussed the possibility that radio pulses 
may be due to hot plasma corotating at discrete positions 
on the light cylinder (e.g. Gold 1968, 1969; Bartel 1978; 
Cordes 1981; Ferguson 1981). These models have since been 
abandoned. Our present model is different in several respects 
(our plasma is cold in its rest frame, the fundamental syn- 
chrotron radiation frequency is the neutron star rotation 
frequency, pulses are due to the 3D structure of the magne- 
tosphere, etc.). 

We would like to conclude this presentation with two 
final points of interest. The first has to do with one of 
the few robust facts about pulsar radio emission, namely 
the characteristic S-curve polarization sweeps. In order for 
our model to survive as a valid alternative to highly suc- 
cessful phenomenological models such as the rotating vec- 
tor model of Radhakrishnan & Cooke (1969), it has to be 
able to explain this fact. This analysis requires knowledge 
of the 3D structure of the pulsar magnetosphere obtained 

2 Our estimates are obtained as follows: for each catalogued pul- 
sar with available radio brightness data we calculate the ON 
brightness as P*/W times the average brightness, where W is 
the observed pulse width. The total radio luminosity is then cal- 
culated in the context of the pulsar synchrotron model in which 
radiation is emitted within ±1/7 radians above and below the 
equator. This is very different from the canonical model where 
radio emission is produced in two beams originating in the polar 
caps. 
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only recently by Spitkovsky (2006) and Kalapotharakos & 
Contopoulos (2009). Let us, however, present here some pre- 
liminary results, and let us introduce the three basic angles 
of the 3D problem, namely the angle a between the rota- 
tion axis and the stellar magnetic moment, the angle i be- 
tween the observer line of sight and the rotation axis, and 
the observed polarization angle ip with tp = along the 
rotation equatorial plane. Let us also introduce the angle 
£ between the rotation equatorial plane and the direction 
connecting the central star with the instantaneous position 
of the tip of the dead zone on the light cylinder. In the con- 
text of our present pulsar synchrotron model, radiation is 
observed whenever the instantaneous direction of the an- 
nular current at the tip of the dead zone happens to lie 
close to our line of sight. The instantaneous direction of 
polarization ip lies along the direction of the instantaneous 
radius of curvature which, we argue, must be close to £. 
Therefore, tp ~ tan _1 [tan(£) sin(i)]. In fig. (4) we show the 
variation of £ with pulsar longitude as obtained in the solu- 
tion of Kalapotharakos & Contopoulos (2009) for inclination 
a = 60°. We also plot the resulting polarization sweep for 
an observer inclination i = 30°. This preliminary result is 
promising enough to justify the continuation of our effort in 
paper II of this series. 

The second point of interest has to do with the equato- 
rial return current sheet beyond the light cylinder. One can 
directly check that this is positively (negatively) charged, 
that is with a surface charge opposite to that of the sepa- 
ratrix current sheet inside the light cylinder. And yet, the 
electric current in both current sheets flows in the same di- 
rection, that is towards (away from) the star. Given the 
limitation on the available number density of charge car- 
riers imposed by the force-free approximation (see discus- 
sion in § 2), in order to supply the inward flowing electrons 
(positrons) and the outward flowing positrons (electrons), 
we need a source of pairs at the tip of the separatrix at the 
position of the light cylinder. This may either be supplied 
as a steady stream of pairs along the separatrix from the 
surface of the star, or through pair creation by some yet to 
be determined process at that position. Note that as these 
particles move away from the tip of the dead zone, they ei- 
ther move along the equator in a trajectory which becomes 
asymptotically radial (i.e. with a radius of curvature much 
larger than r ;c ), or they move inwards where the corotation 
Lorentz factor becomes smaller and smaller. In both cases, 
their contribution to the total synchrotron radiation energy 
content is expected to be much smaller than the radiation 
they emit as they traverse the annular volume of thickness 
d just inside the light cylinder that we discussed in the pre- 
vious section. 
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Figure 1. Schematic of the dependence of the poloidal (dashed 
line) and toroidal (solid line) magnetic field components on radial 
distance R from the star along the equator near the tip of the dead 
zone. The poloidal field varies as oc 1/B? near the star, reaches 
a minimum value at 87% of the light cylinder distance, grows to 
a very large value at the tip of the dead zone as given by eq. \E\ 
and drops to zero outside. The toroidal field is zero inside the 
dead zone, and rises to some finite value outside. The jumps in 
Bp and at the tip of the dead zone take place over some finite 
distance d <C ri c . 



Figure 2. Schematic of electron and positron trajectories along 
the separatrix and the equatorial current sheet respectively for 
aligned magnetic and rotation axes, and detail of the tip of the 
dead zone. The light cylinder is denoted with a dashed line, d ~ 
10 -3 cm is the separatrix thickness. A ~ 1 m is the electron path 
length in the pulsar synchrotron. 
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Figure 3. Pulsar radio luminosity as a function of the pulsar 
period. Open circles: calculated according to eq. 1211 Crosses: es- 
timated from the observed emission at 400 MHz adjusted in the 
context of the pulsar synchrotron model where radiation is emit- 
ted within ±1/7 radians above and below the equator instead of 
in two beams originating in the polar caps. Data are taken from 
the ATNF pulsar catalogue. 



Figure 4. Distribution of £ (dotted line) and estimate of polar- 
ization angle if) (solid line) as functions of the pulsar longitude 
obtained from the numerical solution of Kalapotharakos & Con- 
topoulos (2009) for a = 60° and i = 30°. 



